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The density fluctuation spectrum N(k, u) is calculated for the t-J model in the low-doping regime using a 
slave-fermion method for the constrained fermions. The obtained results for N(k, ui) are in good agreement 
with diagonalization results. The density response is characterized by incoherent, momentum dependent spectral 
functions reaching up to energies ~ 8t and a low-energy structure at energy ~ J due to transitions in the 
quasiparticle band. N(k,u)) is shown to lead to a strong renormalization of planar bond-streching and breathing 
phonon modes with a large phonon linewidth at intermediate momenta caused by the low-energy response. Our 
results are consistent with recent neutron scattering data, showing the peculiar behavior of these modes. 



1. INTRODUCTION 

The density response of doped Mott-Hubbard 
insulators like the high-T c superconductors is ex- 
pected to differ strongly from that of weakly cor- 
related metals. Due to the Mott-Hubbard gap 
the low-energy density response in the doped 
systems is proportional to the number of holes 
and not proportional to the electron concentra- 
tion. Experimentally the peculiarities of the den- 
sity response in high-temperature superconduc- 
tors show up (i) in the anomalous mid-infrared 
absorption in the frequency dependent conductiv- 
ity, (ii) in the renormalization and strong increase 
of linewidth of certain phonon modes upon dop- 
ing, and (iii) the charge instability leading to the 
stripe phase. While static stripes in cuprates are 
found only in the LTT-phase it is expected that 
dynamic stripe like correlations exist also in the 
other cuprate materials. 

Exact diagonalization studies of the dynam- 
ical density response N(k,u>) at large momen- 
tum transfer have revealed several features un- 
expected from the point of view of weakly cor- 
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related fermion systems: (a) the very different 
form of iV(k, lo) compared to the spin response 
function S(k, to), which share common features in 
usual fermionic systems; (b) a strong suppression 
of low energy 2kp scattering in the density re- 
sponse, (c) a broad incoherent peak in the higher 
energy range t < lj < 8t with large dispersion, 
whose shape is rather insensitive to hole concen- 
tration (5 = 0.1 and 0.25) and exchange interac- 
tion J. Furthermore, exact diagonalization stud- 
ies show (d) a low-energy structure on scale J for 
certain momenta which grows in intensity with J. 

While considerable analytical work has been 
done to explain the spin response of the t-J 
modelQ, only few authors analyzed N(k,cj). 
Wang et al. || studied collective excitations in the 
density channel and found sharp peaks at large 
momenta corresponding to free bosons. Similar 
results were obtained by Gehlhoff and Zeyherj| 
using the X-operator formalism. In these calcu- 
lations fluctuations were treated in the leading 
order of a 1/A^-expansion which is however not 
sufficient to explain the incoherent character of 
N(k,uj). It was shown subsequently by Khali- 
ullin and Horsch who used a slave-boson ap- 
proach, that the next order in 1/iV-expansion is 
crucial to obtain renormalized polaron-like boson- 
propagators for the holes, which explained both 
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the incoherent structure at high energy and the 
appearance of a low energy peak associated with 
the coherent polaron motion. The characteristic 
energy scale of the low-energy polaron peak is J ~ 
X J+St, where \ is the uniform RVB-order param- 
eter relevant for the overdoped regime studied in 
that work. In addition on the same energy scale 
J there is a smooth Fermi-liquid like particle-hole 
continuum, which reflects the Fermi surface, with 
total weight ~ <5 2 . Further studies analysed 
the incoherent part of iV(k, ui) yet did not obtain 
the low-energy structure. 

The aim of the present work is to investigate 
(q, ui) in the low-doping regime. We confine 
our study to the AF long-range ordered phase 
which implies a different Fermi surface and col- 
lective spin wave excitations. There are two key 
questions: (i) what theory is needed to produce 
spectra at high energy similar to the large doping 
regime as required by diagonalization studies and 
(ii) what are the implications of the long-range 
magnetic order and the small hole-pocket Fermi 
surface for the low energy charge response? 

2. MODEL AND DENSITY RESPONSE 

We consider here the 2D t — J model as the 
generic model for the description of strongly cor- 
related electrons in cuprates: 



H t j (fe + H.c.) + J ^ SiS s , (1) 

(«',i> 

where the operator Ci a = c, CT (l — n^) describes 
the annihilation of an electron with spin a at site 
i with the constraint of no double occupancy, Si is 
the spin operator, and the summation runs 
over nearest-neighbor bonds on the square lattice. 

We rewrite the constrained electron operators 
in terms of spinless fermions fi which denote the 
holes and use Holstein-Primakoff bosons hi to de- 
scribe the spin degrees of freedom. After transfor- 
mation to momentum representation the Hamil- 
tonian H t j = H t + Hj takes the form 

Ht = ztJ2 tf_ q /k(7k- q b+ + Tk6- q ) + H.c, (2) 



J 



H ' = o Y,^K b \ + + Kki (3) 






Figure 1. Diagrams included in the calculation 
of the memory function. Solid and wavy lines rep- 
resent renormalized spinless fermion and magnon 
propagators, respectively. 



where 7k = 0.5(cosfc x + cosk y ) while z = 4 de- 
notes the number of nearest neighbours. This 
Hamiltonian has the structure of a polaron prob- 
lem where spinless fermions (holes) are coupled 
to collective spin excitations ||. 
The density fluctuation spectrum 



N{k,uj) = -ilmx(k,w 



(4) 



where x(k, uj) = ((n k | n_ k )) w and n k 

E 

relation function 



, a c i a c i a c""-' , is defined via the density cor- 



((n k | n_ k )) u 



[n k (t),n k (0)]>dt. (5) 



Calculation of x(k, ui) by expansion in a small 
parameter usually fails due to the singular char- 
acter of this quantity for small lj. This problem 
can be resolved by rewriting %(k, ui) in terms of 
the static susceptibility Xo(k) and a selfenergy 
M(k,^)[|0]: 

M(k,w) 



X( k ,^) = Xo(k)- 



(6) 



u) + M (k, bj) ' 

The evaluation of the memory function M(k, ui) 
by expanding in a small parameter (here the cou- 
pling constant t) gives correct results for x(k, w) 
valid for all frequencies. From this expression it 
is obvious that an evaluation of M to low order in 
the coupling constant t implies for the charge sus- 
ceptibility x a partial resummation up to infinite 
order. We shall find that this form is crucial for 
the proper description of the charge excitations 
at high energy, which are of collective nature and 
induced by the denominator uj+M which can be- 
come small. In second order approximation in t 
we can write: 



wx(k,w) ~xo(k)M(k,w), 



(7) 
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Figure 2. Density fluctuation spectra for 6 = 0.1 
calculated via the slave-fermion approach for J = 
0.4 and momenta along (tt,tt) and (7r,0) direc- 
tions. The energy unit is t = 1. 



and determine the memory function by a low or- 
der calculation of x- Repeated application of 
equation of motions yields 



M(k,w) = [/(k,w)-/(k,0)]/xo(k), 



(8) 



where /(k,w) = — ((rik | ^-k))w The diagram- 
matic representation of the /(k, w) function is 
shown in Fig. 1. The fermion propagators are 
evaluated within selfconsistcnt Born approxima- 
tion at finite doping concentration; more details 
will be given elsewhere. The still unknown xo(k) 
is determined from the integral relation 



iV(k) = / N(k,u)dbJ, 



(9) 



and the fact that the static charge structure factor 
N(k) is well approximated in the strong correla- 
tion regime by the N(k) for spinless fermions|]lT|. 

The density fluctuation spectra obtained by the 
slave-fermion calculation are displayed in Fig. 2. 
N(k, uj) is characterized by a broad continuum 
which moves to higher energy as k increases. At 
(7r, 7r) the spectrum peaks at to ~ Qt, i.e., the 
variation of iV(k, u>) with k and the energy scale 



are fully consistent with the diagonalization re- 
sults for small systems 0. Low energy structures 
on scale J, due to quasiparticle scattering be- 
tween hole pockets, appear at small and interme- 
diate momentum transfers and vanish at the zone 
boundary. The pronounced low-energy structure 
at large momentum transfer in (n, 0) direction 
found in the spin liquid regime || is absent. 

3. BOND-STRETCHING PHONONS 

The high energy planar breathing and bond- 
stretching oxygen modes in (ir, ir) and (ir, 0) di- 
rection, respectively, show the strongest dop- 
ing dependence from all phonon modes |l2| 
which appears to be a generic feature of cuprates. 
Inelastic neutron scattering studies of the LO 
mode in (7r,0) direction in L^-zSr^CuC^ (x — 
0.15)|l2|,|l3| reveal a ~ 10% renormalization near 
(tt, 0) as compared to the undoped system, while 
there is no significant renormalization at the T- 
point. Surprisingly the lincwidth assumes a max- 
imum at intermediate k|l3[]. 

The coupling of these phonon modes to density 
fluctuations is described by pi 



He-ph = gJ2 fi 



(10) 



where u l x , u l y are the displacements of the four 
oxygen neighbors of the hole, g is the electron- 
phonon coupling constant. This coupling origi- 
nates from the change of the Zhang-Rice singlet 
energy by breathing and bond-stretching modes. 

The phonon Green's function has the following 
form for k = (k, 0) 



D(k,uj) 



u> 



+4£sin^fx(fc,u;)) 



(11) 



here we assume a k-independent bare phonon 
frequency loq = 0.3 and a normalized coupling 
constant £ = g 2 /Kt = 0.25 |lj]. For momenta 
k = (k,k) the constant 4£ in Eq.(ll) should be 
replaced by 8£. The phonon spectral function 
B(k,uj) = —^ImD(k,uj) is displayed in Fig. 3. 
Figure 3 shows a large increase of the renormal- 
ization of Wk with increasing k, which is similar 
for the two directions, with rather abrupt changes 
in cjk at intermediate k. The damping appears to 
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Figure 3. Phonon spectral function at 6 = 0.1 
for different momenta along (k, 0) and (k, k) di- 
rections, where k = ir/8 for the top panel and 
increases in steps of 7r/8 in downward direction. 
Here loq = 0.3; arrows indicate u>k along [k, k). 



be largest not at the zone-boundary but at in- 
termediate k. We note that earlier calculations 
for larger doping, S = 0.15 — 0.25, show a sig- 
nificantly larger energy change at (jr, 0) than at 
(-7T, 7r), which was due to the pronounced low en- 
ergy polaron related structure in N(\c,u>) near 
(ir, 0) in the spin disordered stateQ. The low- 
energy structure of N(k,uj) in the AF ordered 
system considered here is markedly different, and 
suggests nontrivial changes for the phonon renor- 
malization as function of doping. 

Hence we conclude by stressing that phonon 



modes, and in particular those discussed here, 
provide a sensitive probe for the anomalous low- 
energy charge response of doped Mott insulators. 
Furthermore it is expected that measurements at 
different doping concentrations will reveal signif- 
icant changes of phonon energies and damping 
due to the changes in the spin response and in 
the topology of the Fermi surface, which in turn 
strongly affect the charge response at low energy. 
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